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Abstract:  In  2003,  the  U.S.  Army  Corps  of  Engineers  issued  ERDC/ITL 
TR-03-3,  which  interpreted  eight  of  the  44  full-scale,  low-velocity, 
controlled-impact  barge  train  experiments  conducted  at  the  decom¬ 
missioned  Gallipolis  Lock  at  Robert  C.  Byrd  Lock  and  Dam.  An  easy-to-use 
“empirical  correlation”  was  derived,  reporting  the  maximum  impact  force 
(normal  to  the  wall)  as  a  function  of  the  linear  momentum  normal  to  the 
wall  (immediately  before  impact).  This  empirical  correlation  was  cited  in 
ETL 1110-2-563  (April  2004)  for  impacts  with  stiff-to-rigid  walls  that  do 
not  involve  damage  to  either  the  corner  barge  or  the  wall;  a  limiting 
impact  force  of  800  kips  (capping  the  empirical  correlation)  was  cited.  In 
this  technical  report,  we  provide  a  basis  for  revising  the  “capping  force” 
cited  in  ETL  1110-2-563,  based  on  the  computation  of  the  “crushing”  force 
imparted  during  a  glancing  blow  of  the  impact  corner  at  the  bow  with  a 
lock  approach  wall. 

A  second  limiting  force  due  to  the  limit  state  of  lashing  forces  during  a 
glancing  blow  impact  is  addressed  in  the  USACE  technical  report 
ERDC/ITL  TR-05-1.  The  lower  of  the  two  limiting  forces— due  to  the 
crushing  of  the  barge  impact  corner  or  to  the  failure  of  the  lashings— 
“caps”  the  empirical  correlation. 
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1  Introduction 

1.1  Background:  barge  train  maximum  impact  forces 

Locks  are  a  necessary  structural  feature  found  at  every  dam  within  the 
U.S.  inland  waterways  navigation  system.  This  network  of  rivers  is  an 
essential  component  of  the  nation’s  transportation  infrastructure  system,  a 
system  key  to  national  commerce.  Locks  allow  for  groups  of  barges,  lashed 
together  to  form  barge  trains,  to  negotiate  the  changes  in  river  elevation  at 
the  dams.  One  of  the  most  frequent  loads  applied  to  the  locks  of  the 
U.S.  inland  waterway  system  is  the  impact  made  by  a  barge  train  on  the 
approach  walls  as  the  barge  train  aligns  itself  to  transit  the  lock.  Conse¬ 
quently,  this  load  case  represents  one  of  the  primary  design  loads  con¬ 
sidered  for  lock  approach  walls. 

This  research  report  discusses  the  results  of  a  series  of  nonlinear  finite  ele¬ 
ment  analyses  computing  the  limiting  impact  force  due  to  yielding  and 
buckling  of  the  deck  and  skin  plates  and  the  internal  structural  frame  within 
the  “impact  corner”  of  the  bow  of  a  jumbo  open-hopper  barge  during  a 
glancing  blow  impact  with  a  lock  approach  wall.  The  “structural  concept” 
can  be  explained  as  follows.  Because  of  the  elasto-plastic  and  limiting  strain 
material  characteristics  of  steel— combined  with  the  structural  layout  of  the 
deck  and  hull  plate,  the  internal  structural  plates,  and  the  angle  steel  of  the 
internal  trusses— the  multi-degree-of-freedom  structure  of  the  barge  bow 
provides  a  limiting  force  resistance  during  an  impact  event.  When  placed  in 
a  severe-impact  environment,  in  which  the  impact  corner  of  the  bow  begins 
to  “crush,”  the  barge  bow  will  act  like  a  structural  “fuse-plug”  to  provide  for 
a  limiting  impact  force  applied  to  the  approach  wall  by  the  barge  train. 

Thus,  there  is  an  upper-bound  force  that  approach  walls  will  be  subjected  to 
during  a  glancing  blow  impact  event.  This  report  summarizes  the  research 
effort  investigating  the  magnitude  of  this  limiting  force.  The  barge  used  in 
this  study  is  a  rake  section,  jumbo  open-hopper  barge  200-ft  long,  35-ft 
wide,  and  13-ft  high  at  the  hopper  section  region. 

1.1.1  Empirical  correlation 

Lock  approach  walls  are  designed  for  usual,  unusual,  and  extreme  loads.  In 
the  past,  the  primary  focus  of  engineers  performing  impact  computations 
has  been  on  the  lock  approaches  where  the  worst-case  events  and  damage  to 
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the  barge  train  and/or  wall  are  likely  to  occur  (e.g.,  extreme  loads  where 
vessel  control  is  lost).  The  U.S.  Army  Corps  of  Engineer’s  (USACE)  initial 
guidance,  engineer  technical  letter  ETL 1110-2-338,  for  computing  glancing- 
blow  maximum  impact  forces  for  barge  trains  impacting  approach  walls  was 
rescinded  in  2001.  The  impact  force  computations  in  ETL  1110-2-338  were 
based  on  crushing  of  the  impact  corner  of  the  barge,  a  response  that  is  usu¬ 
ally  associated  with  an  extreme  load  case,  and  not  the  no-damage  usual  load 
case,  nor  the  minor  damage  unusual  load  case.  In  April  2004,  ETL  1110-2- 
563  “Barge  Impact  Analysis  for  Rigid  Walls”  was  issued  by  Headquarters, 
USACE.  The  impact  force  computations  are  based  on  the  “empirical  correla¬ 
tion”  developed  by  Arroyo  et  al.  (2003),  shown  in  Figure  1.1.  Details  regard¬ 
ing  its  development  are  summarized  in  Appendix  A. 


0  200  400  600  800  1000  1200 

Linear  Momentum  Normal  to  the  Wall  (k-s) 

(mVnorm) 

Figure  1.1.  Empirical  correlation  using  the  linear  momentum  normal  to  the  wall  concept; 
from  Figure  6.3  in  Arroyo  et  al.  (2003). 

This  figure  reflects  advances  in  the  computations  of  usual  and  some  unusual 
impact  design  loads  during  which  no  damage  occurs  to  the  wall  and  the 
barges,  and  no  failure  occurs  in  the  lashings  that  bind  the  barges.  The 
reduction  of  the  Patev  et  al.  (2003)  load  data  measured  during  full-scale, 
low-velocity,  controlled-impact  barge  experiments  using  a  3-by-5  (three 
wide  and  five  long)  15-barge  train  impacting  a  wall  at  a  decommissioned 
USACE  lock  is  discussed  in  Arroyo  et  al.  (2003),  along  with  the  resulting 
empirical  correlation.  Figure  1.1  idealizes  the  interrelationship  between 
maximum  impact  force  imparted  to  the  wall  and  the  linear  momentum 
(immediately  before  impact)  of  the  barge  train  normal  to  the  wall.  Linear 
momentum  conveniently  represents  two  key  demand  variables:  the  mass  of 
the  barge  train  and  its  velocity  normal  to  the  wall. 
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The  empirical  correlation  in  Figure  1.1  is  shown  to  be  unbounded.  In  actu¬ 
ality,  it  is  bounded  by  a  limiting  force  that  functions  as  an  asymptote  to  the 
empirical  correlation,  as  idealized  in  Figure  1.2. 


Figure  1.2.  Asymptote  to  the  empirical  correlation. 


This  limiting  impact  force  in  Figure  1.2  results  from  either: 

•  the  failure  of  the  lashings  that  bind  the  barges  of  the  barge  train 
together,  or 

•  the  plastic  yielding  of  the  structural  members  (i.e.,  yielding  and  buckling 
of  the  deck  and  skin  plates  and  the  internal  structural  frame)  in  the  cor¬ 
ner  of  the  barge  that  impacts  the  approach  wall. 

Whichever  limiting  impact  force  value  is  less  becomes  the  limiting  force  to 
the  empirical  correlation  in  Figure  1.1.  A  limiting  impact  force  of  800  kips, 
based  on  engineering  judgment,  was  cited  in  the  2004  initial  version  of 
ETL 1110-2-563. 

Since  2004,  there  has  been  additional  research  into  the  limiting  impact 
force  due  to  the  failure  of  the  lashings  for  various  sizes  of  barge  trains. 
Arroyo  and  Ebeling  (2004)  developed  and  fully  described  the  theory  for  two 
probable  idealized  lashing  failure  mechanisms  during  a  glancing  blow 
impact  with  an  approach  wall.  They  are  described  as  the  transverse  and 
corner  failure  mechanisms.  This  formulation  was  implemented  by  Arroyo 
and  Ebeling  in  the  PC  software  program  Limit  LASHING.  Arroyo  and 
Ebeling  (2005)  present  a  full  assessment  and  a  parametric  study  of  these 
failure  mechanisms.  The  parametric  study  is  used  to  quantify  the  magnitude 
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of  the  maximum  impact  force  as  a  function  of  a  range  of  values  of  each  of 
the  primary  variables  within  the  barge  train  system.  This  parametric  study 
was  based  on  a  numerical  model  for  computing  barge  impact  forces  based 
on  the  ultimate  strength  of  the  lashings  between  barges,  computed  using  the 
personal  computer  software  program  Limit_LASHING.  This  numerical 
study  resulted  in  the  computation  of  limiting  impact  forces  below  and  above 
the  8oo-kips  limiting  force  cited  in  the  2004  version  of  ETL 1110-2-563. 

The  value  for  the  limiting  force  due  to  lashing  failures  was  dependent  upon 
the  size  of  the  barge  train,  the  size  of  the  wire  rope  (i.e.,  lashing),  the  orien¬ 
tation  of  the  lashing  layout  connecting  the  barges,  the  number  of  wraps 
around  the  bits,  and  whether  the  rope  was  new  or  used.  The  results  of  the 
numerous  parametric  analyses  are  presented  in  graphical  form  in  Arroyo 
and  Ebeling  (2005). 

This  research  report  extends  the  limiting  impact  force  due  to  lashing  fail¬ 
ures  to  a  discussion  of  results  for  a  series  of  LS-DYNA  nonlinear  finite  ele¬ 
ment  analyses,  computing  the  limiting  impact  force  due  to  yielding  and 
buckling  of  the  deck  and  skin  plates  and  the  internal  structural  frame  in  the 
front  corner  (e.g.,  Figure  1.3)  of  a  jumbo  open-hopper  barge  that  impacts 
the  wall.  With  the  results  of  the  Arroyo  and  Ebeling  (2005)  parametric 
study  and  those  contained  in  this  report,  the  limiting  impact  force  for  the 
empirical  correlation  may  now  be  established  for  a  barge  train  of  a  given 
size. 

1.1.2  Yielding  and  plate  buckling  analyses 

In  the  LS-DYNA  nonlinear  finite  element  analyses  discussed  in  this  report, 
the  impact  corner  of  the  barge  and  the  approach  wall  (modeled  as  a  non¬ 
penetrating  “rigid”  structure)  are  brought  into  contact  with  one  another  at  a 
constant  velocity  in  the  numerical  simulation  of  the  “crushing”  of  the 
impact  corner  of  the  barge.  “Crushing”  of  the  elliptical  impact  corner  con¬ 
tinues  until  36  in.  of  penetration  of  the  approach  wall  into  the  barge  impact 
corner  is  achieved.  Figure  1.4  (top)  shows  an  overhead  view  of  the  barge 
impact  corner  at  initial  contact  with  the  wall  in  the  numerical  analysis  using 
LS-DYNA;  Figure  1.4  (bottom)  shows  an  overhead  view  of  the  barge  impact 
corner  after  36  in.  of  penetration  for  one  of  the  numerical  evaluations. 
Figure  1.5  shows  the  corresponding  resultant  contact  force  normal  to  the 
“rigid”  wall,  as  computed  by  LS-DYNA  during  the  course  of  this  numerical 
analysis.  The  key  observation  is  that  the  contact  force 
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Figure  1.3.  View  of  the  impact  corner  of  the  jumbo 
open-hopper  barge  used  in  this  numerical  study. 


between  the  “rigid”  nonpenetrating  approach  wall  and  the  elliptical  impact 
corner  increases  in  magnitude  with  wall  deformation  into  the  barge  up  to  a 
peak  force  value:  the  limiting  contact  force.  This  limiting  contact  force  value 
becomes  the  asymptote  to  the  empirical  correlation  in  Figure  1.1,  as  ideal¬ 
ized  in  Figure  1.2. 

Figure  1.5  is  an  example  of  typical  resultant  contact  force  normal  to  the 
“rigid”  (approach)  wall  versus  permanent  deformation  computed  at  the 
impact  corner  of  the  barge  in  the  LS-DYNA  nonlinear  numerical  analyses. 

In  this  figure,  the  normal  force  ranges  in  value  from  o  kips  to  a  maximum 
force  of  3,497  kips  at  4.8  in.  of  permanent  normal  deformation  at  the 
impact  corner.  All  LS-DYNA  analyses  were  continued  out  to  36  in.  of  per¬ 
manent  deformation  of  the  rigid  wall  into  the  impact  corner  of  the  barge. 
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Top  Edge  of 
Hopper  Area 


pjgr. 


Figure  1.4.  Overhead  views  of  the  barge  impact  corner  at  initial  contact  with  the  wall  (top) 
and  after  36  in.  of  penetration  (bottom). 
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Figure  1.5.  Force  versus  displacement  of  the  impact  corner  for  Case  No.  9.  Unusual  loading 
condition  with  an  approach  angle  of  10°. 

1.2  Nonlinear  structural  dynamic  analyses  using  LS-DYNA 

LS-DYNA  is  a  general  purpose  transient  finite  element  program  that  is  used 
for  analyzing  complex  structural  dynamics  problems,  such  as  the  crushing 
of  the  impact  corner  of  the  bow  of  a  barge  during  a  glancing  blow  impact 
with  an  approach  wall.  It  uses  a  central  difference  scheme  to  solve  the  equa¬ 
tion  of  motion  in  time.  Key  features  used  in  the  analyses  discussed  in  this 
report  are  nonlinear  dynamics,  use  of  a  nonlinear  elastic-plastic  constitutive 
model  for  all  (A-36)  steel  plates  and  angles  comprising  the  barge,  and  use  of 
contact  surface  formulation  for  a  flexible-body-to-rigid-body  contact. 
Lagrangian  shell  elements  are  used  to  model  the  barge  components. 
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1.3  Report  contents 

Chapter  2  discusses  the  finite  element  mesh  for  the  bow  of  the  jumbo  open- 
hopper  barge  assembled  using  TrueGRID  (Rainsberger  2006)  and  modeled 
in  the  nonlinear  LS-DYNA  analyses  as  well  as  the  assigned  material  proper¬ 
ties.  Material  properties  corresponding  to  A-36  steel,  obtained  from  tests 
conducted  on  plate  steel  from  a  barge,  are  assigned  in  these  numerical 
analyses. 

Chapter  3  reports  on  the  results  of  the  24  LS-DYNA  nonlinear  finite  element 
analyses.  The  approach  angle  of  the  barge  with  the  approach  wall,  as  well  as 
the  approach  velocity,  is  varied  among  the  analyses. 

Summary  and  conclusions  are  contained  in  Chapter  4. 

Appendix  A  discusses  the  background  for  the  full-scale,  low-velocity 
controlled-impact  barge  experiments  and  the  development  of  the  empirical 
correlation. 

Appendix  B  discusses  the  results  of  an  evaluation  of  influence  of  the  number 
of  processors  and  the  accuracy  of  computed  results  for  the  nonlinear  analy¬ 
ses  using  fully  integrated  shell  elements  and  reduced  integration  elements 
with  hourglass  control. 
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2  Finite  Element  Mesh  and 
Material  Properties 

2.1  Introduction 

This  chapter  discusses  the  development  of  the  finite  element  mesh  and  the 
material  properties  assigned  to  the  structural  steel  comprising  the  barge. 
The  object  of  the  nonlinear  finite  element  analysis  (using  LS-DYNA)  is  to 
compute  the  limiting  impact  force  due  to  the  yielding  and  buckling  of  the 
plates  and  internal  structural  framing  at  the  impact  corner  of  the  barge 
during  its  glancing  blow  impact  with  a  lock  approach  wall;  therefore,  only 
the  bow  region  of  the  jumbo  open-hopper  barge  needs  to  be  modeled  in  the 
numerical  analysis  (Figure  2.1,  raked  section).  Additionally,  because  the 
nonlinear  deformations  are  concentrated  in  the  front  half  of  the  bow  region 
closest  to  the  impact  corner,  only  half  the  bow  is  modeled  in  the  analysis. 


Figure  2.1.  Elliptical  impact  corner  and  bow  of  a  raked  section,  jumbo  open-hopper  bow. 
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This  half-bow  section  is  17.5-ft  wide  and  27.7-ft  long.  The  bow  is  15.5-ft 
high;  at  the  intersection  of  the  inclined  open-hopper  face  with  the  flat 
hopper  base,  the  barge  is  13-ft  high. 

2.2  Finite  element  mesh 

The  finite  element  software  TrueGRID  is  used  to  construct  the  finite  ele¬ 
ment  mesh  of  the  bow  of  the  jumbo  open-hopper  barge.  Fully  integrated 
LS-DYNA  shell  elements  (*SECTION_SHELL;  ELFORM  16)  are  used  in  a 
detailed  model  (Figure  2.1)  of  the  jumbo  open-hopper  barge  bow  structural 
steel  deck,  hull,  and  internal  structural  plates,  as  well  as  the  internal  angle 
steel  comprising  the  internal  structural  trusses. 1  The  mesh  consists  of 
179,238  nodes  and  177,241  elements.  The  resulting  mesh  of  the  front  half  of 
the  bow  of  the  barge  is  shown  in  Figures  2. 2-2. 5.  Also  identified  in  these 
figures  (by  color  designation)  are  the  thicknesses  of  the  structural  deck, 
hull,  and  hopper  plates  of  the  bow  of  the  jumbo  open-hopper  barge.  The 
headlog  and  elliptical  impact  corner  sections  of  the  bow  are  shown  in  yellow 
in  Figures  2.3  and  2.5.  The  rounded  hull  plate  connecting  the  front  and  side 
hull  plates  (in  blue)  is  shown  in  brown  in  Figures  2.2  and  2.3.  With  the 
mesh  terminating  along  the  center  line  of  the  bow  of  the  barge,  Figure  2.4 
allows  for  a  view  of  the  internal  truss  system  within  the  bow  of  the  barge  as 
well  as  the  deck  and  hull  plates.  There  are  seven  internal  trusses  (six  “stan¬ 
dard”  longitudinal  and  one  outboard  longitudinal  truss)  contained  within 
this  half-bow  structural  section  model.  The  trusses  are  made  out  of  A-36 
steel  angles.  Figure  2.6(a)  provides  an  internal  view  of  the  actual  structural 
trusses  and  elliptical  (impact)  corner  at  the  bow  of  the  barge  whereas  Fig¬ 
ure  2.6(b)  shows  the  corresponding  view  of  the  structural  model.  Figure  2.7 
(left)  shows  the  view  of  the  elliptical  corner  during  its  assemblage  and  with¬ 
out  the  deck  plate,  whereas  Figure  2.7  (right)  shows  the  corresponding 
mesh.  Figure  2.8  shows  the  internal  structural  members  with  the  deck  and 
hull  plates  removed. 


1  LS-DYNA  reduced  integrated  shell  elements  (*SECTION_SHELL;  ELFORM  =  2;  NUMBER  OF  INTEGRATION 
POINTS  =  3)  with  hourglass  control  (*FIOURGLASS;  Flourglass  Control  Type  =  4,  Flanagan-Belytschko 
stiffness  form)  were  tried  to  save  on  execution  time.  Reduced  integration  resulted  in  execution  times 
nearly  one-third  that  for  the  fully  integrated  elements.  Flowever,  the  results  obtained  varied  significantly  if 
the  number  of  processors  used  to  run  a  problem  was  varied.  Therefore,  the  fully  integrated  shell  element 
model,  used  as  the  result  obtained  from  runs  with  varying  numbers  of  processors,  was  observed  to  be 
more  consistent. 
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Figure  2.2.  Mesh  of  bow— side  view  of  the  hull  and  hopper  plates. 


Figure  2.3.  Front  elliptical  corner  view  of  the  deck  plates,  headlog,  and  elliptical  corner  of  the  bow. 
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Figure  2.4.  View  of  the  bow  of  the  barge  at  the  centerline  cut  of  the  bow. 


Figure  2.5.  Overhead  view  of  the  bow  of  the  barge. 
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Figure  2.6.  Internal  views  of  the  structural  trusses  and  elliptical  corner  at  the  bow  of  the  barge. 
(Top)  Actual  internal  structural  truss  members  and  the  elliptical  corner.  (Bottom)  Corresponding 
finite  element  model  and  material  zone  steel  members,  identified  by  steel  thicknesses. 
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Figure  2.7.  Views  of  the  elliptical  corner  structural  member.  Interior  view  of  the  assemblage  of 
the  elliptical  corner  (left)  resting  on  the  ground  from  a  bow  point  of  view.  (Right)  Corresponding 

finite  element  model  view. 
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Figure  2.8.  Internal  structural  members  with  the  deck  and  hull  plates  removed. 

(Top)  Centerline/elevated  view  of  internal  members  (elliptical  corner  plates  in  green  on  right). 
(Bottom)  Centerline/elevated  view  of  internal  members  with  three-quarter-inch-thick  plate  on 

elliptical  corner  and  headlog. 
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2.3  Material  modeling  and  material  properties 

With  regard  to  material  modeling  used  for  the  structural  steel,  significant 
yielding  is  expected  to  occur  during  crushing  of  the  deck,  hull  plates,  and 
the  internal  structural  members  in  the  impact  corner  zone  of  the  bow.  As 
such,  all  shell  elements  in  the  model  of  the  bow  were  specified  with  a 
plastic-multilinear  material  model  matching  the  stress-strain  data  derived 
from  tests  conducted  on  plate  steel  obtained  for  an  actual  barge.  A  large- 
displacement,  large  strain  formulation  was  used  for  the  shell  elements  in 
the  LS-DYNA  impact  analyses.  Therefore,  the  true  (Cauchy)  stress  and  true 
(logarithmic)  strain  data  was  used  (Figure  2.9).  This  material  model  for 
A-36  steel  was  obtained  from  tests  conducted  at  the  University  of  Florida  on 
standard  18-in.  tension  coupon  (Consolazio  et  al.  2002).  The  steel  model  for 
A-36  steel  has  an  initial  yield  of  36  kips/in.2  (2.48c  +  08  N/m2),  a  yield 
strain  of  0.0012,  and  a  failure  strain  of  0.2.  In  the  LS-DYNA  analyses,  the 
material  model  was  specified  as  an  effective  stress  versus  effective  plastic 
strain,  LS-DYNA  material  number  24, 

*MAT_PIECEWISE_LINEAR_PLASTICITY.  Other  material  model 
parameters  specified  as  input  are: 

•  E  =  30,000  ksi 

•  Poisson’s  ratio  =  0.33 

•  Yield  stress  =  48  ksi 

•  Rupture  stress  =  86.6  ksi  at  25%  strain 

•  20%  plastic  strain  failure  (i.e.,  EPPF). 

Note  that  the  rupture  strain  of  20%  is  specified  in  these  LS-DYNA  impact 
analyses.  The  implications  are  that  when  the  strain  reaches  20%,  the  strain 
integration  point  is  removed  from  the  analysis  by  LS-DYNA. 

The  LS-DYNA  material  number  24  model  allows  for  the  addition  of  strain 
rate  effects  on  the  material  using  the  Cowper  Symonds  model.  However, 
these  simulations  did  not  implement  any  strain  rate  effects  for  the  A-36 
steel  material  because  the  barge  impact  events  are  not  high-speed  impact 
simulations. 
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Figure  2.9.  True  stress  versus  true  strain  for  A-36  structural  steel 
(after  Consolazio  et  al.  2002). 
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3  Limiting  Impact  Force  Computations 

3.1  Introduction 

This  chapter  discusses  the  results  of  a  series  of  24  nonlinear  finite  element 
analyses  using  the  177,241  element  mesh  with  the  assigned  A-36  steel 
material  properties  discussed  in  Chapter  2.  The  object  of  the  nonlinear 
finite  element  analysis  (using  LS-DYNA)  is  to  compute  the  limiting  impact 
force  due  to  yielding  and  buckling  of  the  plates  and  internal  structural 
framing  at  the  impact  corner  of  the  barge  during  its  glancing  blow  impact 
with  a  lock  approach  wall.  Because  the  nonlinear  deformations  are  con¬ 
centrated  in  the  front  half  of  the  bow  region  closest  to  the  impact  corner, 
only  half  of  the  jumbo  open-hopper  bow  shown  in  Figure  2.1  is  modeled  in 
the  analysis. 

3.2  Nonlinear  impact  analyses 

Lock  approach  walls  are  designed  for  usual,  unusual,  and  extreme  loads. 
The  report  “Barge  Impact  Analysis  for  Rigid  Walls”  (ETL 1110-2-563)  pro¬ 
vides  a  summary  of  the  design  requirements  for  these  three  load  cases. 

The  load  cases  are  defined  in  terms  of  the  annual  probability  of  exceed¬ 
ance.  The  design  information  contained  in  the  2004  ETL  report  is  sum¬ 
marized  in  Table  3.1.  This  information  also  includes  ranges  for  nonsite- 
specific  velocities  (expressed  in  local  barge  coordinates)  for  the  three  load 
cases  to  use  in  preliminary  analyses  when  site-specific  traffic  study  results 
are  not  yet  available. 


Table  3.1.  Three  design  load  condition  categories,  frequency  of  loadings,  and  typical  ranges 
for  non-site-specific  impact  angles  and  approach  velocities  (from  ETL  1110-2-563). 


Load 

Condition 

Annual  Probability 
of  Exceedance 
(Return  Period) 

Performance 

Criteria 

Forward 
Velocity,  Vx 
(ft/sec) 

Lateral 
Velocity,  Vy 
(fl/sec) 

Approach 
Angle  0(°) 

Usual 

>0.1  (1-10  years) 

No  damage 

0. 5-2.0 

0.01-0.1 

5-10 

Unusual 

<0.1  but 
>0.00333 
(10-300  years) 

Reparable 

damage 

3. 0-4.0 

0.4-0. 5 

10-20 

Extreme 

<0.00333 
(>300  years) 

Noncollapse 

4. 0-6.0 

<1.0 

20-35 
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The  velocity  of  a  barge  is  usually  specified  in  local  barge  coordinates  as 
shown  in  Figure  3.1  and  is  a  vector  quantity.  The  velocity  of  a  barge  train  is 
usually  specified  in  the  local  barge  axis;  longitudinal  (local  “x”  axis)  and 
transverse  (local  “y”  axis).  In  this  case  two  velocities  are  specified,  Vx  and 
Vy.  To  obtain  the  velocity  normal  to  the  wall,  an  axis  transformation 
equation  is  needed.  This  expression  is  the  following: 


par 


\V„ 


=  c 


-1 


\VX 

\v„ 


(3.1) 


where  [Cl  1 


.  and  Vpar  and  Vnorm  are  the  velocity  parallel 


cos0  -sin0 
sin0  cos0 

(global  “ X ”  axis)  and  normal  (global  “F’  axis)  to  the  wall,  respectively. 
Equation  3.1  can  be  obtained  easily  from  Figure  3.1. 


Using  the  minimum/maximum  ranges  in  Table  3.1  for  non-site-specific 
(1)  impact  (i.e.,  approach)  angles,  (2)  forward  velocities,  and  (3)  lateral 
velocities  for  each  of  the  three  load  cases,  a  combination  of  24  sets  of 
approach  angles  and  velocities  normal  to  the  wall  were  derived.  The 
results  of  these  computations  are  summarized  in  Table  3.2,  along  with 
other  data  pertinent  to  the  LS-DYNA  analyses.  This  information  served  as 
the  basis  for  the  24  LS-DYNA  runs;  run  numbers  are  also  included  in  this 
table. 


Table  3.2.  Twenty-four  LS-DYNA  runs. 


Category 

Impact  Angle 

0(°) 

Forward  Velocity,  Vx 

Lateral  Velocity,  Vy 

Approach  Velocity,  V 

Normal  Velocity 

Parallel  Velocity 

Case 

No. 

(in./sec) 

(ft/sec) 

(in./sec) 

(ft/sec) 

(in./sec) 

(ft/sec) 

P(°) 

(in. /sec) 

(ft/sec) 

(in./sec) 

(ft/sec) 

Usual 

5 

6 

0.5 

0.12 

0.01 

6.0012 

0.5001 

1.1458 

0.6425 

0.0535 

5.9667 

0.4972 

1 

Usual 

5 

6 

0.5 

1.2 

0.1 

6.1188 

0.5099 

11.3099 

1.7184 

0.1432 

5.8726 

0.4894 

2 

Usual 

5 

24 

2 

0.12 

0.01 

24.0003 

2. 

0.2865 

2.2113 

0.1843 

23.8982 

1.9915 

3 

Usual 

5 

24 

2 

1.2 

0.1 

24.03 

2.0025 

2.8624 

3.2872 

0.2739 

23.8041 

1.9837 

4 

Usual 

10 

6 

0.5 

0.12 

0.01 

6.0012 

0.5001 

1.1458 

1.1601 

0.0967 

5.888 

0.4907 

5 

Usual 

10 

6 

0.5 

1.2 

0.1 

6.1188 

0.5099 

11.3099 

2.2237 

0.1853 

5.7005 

0.475 

6 

Usual 

10 

24 

2 

0.12 

0.01 

24.0003 

2. 

0.2865 

4.2857 

0.3571 

23.6145 

1.9679 

7 

Usual 

10 

24 

2 

1.2 

0.1 

24.03 

2.0025 

2.8624 

5.3493 

0.4458 

23.427 

1.9523 

8 

Unusual 

10 

36 

3 

4.8 

0.4 

36.3186 

3.0265 

7.5946 

10.9784 

0.9149 

34.6196 

2.885 

9 

Unusual 

10 

36 

3 

6 

0.5 

36.4966 

3.0414 

9.4623 

12.1602 

1.0133 

34.4112 

2.8676 

10 

Unusual 

10 

48 

4 

4.8 

0.4 

48.2394 

4.02 

5.7106 

13.0622 

1.0885 

46.4373 

3.8698 

11 

Unusual 

10 

48 

4 

6 

0.5 

48.3735 

4.0311 

7.125 

14.244 

1.187 

46.2289 

3.8524 

12 

Unusual 

20 

36 

3 

4.8 

0.4 

36.3186 

3.0265 

7.5946 

16.8232 

1.4019 

32.1872 

2.6823 

13 

Unusual 

20 

36 

3 

6 

0.5 

36.4966 

3.0414 

9.4623 

17.9509 

1.4959 

31.7768 

2.6481 

14 

Unusual 

20 

48 

4 

4.8 

0.4 

48.2394 

4.02 

5.7106 

20.9275 

1.744 

43.4635 

3.622 

15 

Unusual 

20 

48 

4 

6 

0.5 

48.3735 

4.0311 

7.125 

22.0551 

1.8379 

43.0531 

3.5878 

16 

Extreme 

20 

48 

4 

7.2 

0.6 

48.537 

4.0447 

8.5308 

23.1828 

1.9319 

42.6427 

3.5536 

17 

Extreme 

20 

48 

4 

10.8 

0.9 

49.2 

4.1 

12.6804 

26.5656 

2.2138 

41.4114 

3.451 

18 

Extreme 

20 

72 

6 

7.2 

0.6 

72.3591 

6.0299 

5.7106 

31.3912 

2.6159 

65.1953 

5.4329 

19 

Extreme 

20 

72 

6 

10.8 

0.9 

72.8055 

6.0671 

8.5308 

34.7741 

2.8978 

63.9641 

5.3303 

20 

Extreme 

35 

48 

4 

7.2 

0.6 

48.537 

4.0447 

8.5308 

33.4296 

2.7858 

35.1895 

2.9325 

21 

Extreme 

35 

48 

4 

10.8 

0.9 

49.2 

4.1 

12.6804 

36.3785 

3.0315 

33.1247 

2.7604 

22 

Extreme 

35 

72 

6 

7.2 

0.6 

72.3591 

6.0299 

5.7106 

47.1954 

3.9329 

54.8492 

4.5708 

23 

Extreme 

35 

72 

6 

10.8 

0.9 

72.8055 

6.0671 

8.5308 

50.1443 

4.1787 

52.7843 

4.3987 

24 

ro 

o 


ERDC/ITL  TR-08-2 


ERDC/ITL  TR-08-2 


21 


There  are  four  approach  (i.e.,  impact)  angles  used  in  the  24  nonlinear 
finite  element  analyses:  50,  io°,  20°,  and  350.  The  velocity  normal  to  the 
approach  wall  varies  from  a  low  of  0.64  in./sec  (0.054  ft/sec)  to  a  high 
value  of  50.1  in./sec  (4.18  ft/sec).  The  authors  of  this  report  consider  these 
barge  impact  problems  to  be  low-velocity  impacts,  compared  to  “crash” 
analyses  involving  the  impacts  between  vehicles  and  highway  barriers.1 
Each  of  the  24  LS-DYNA  analyses  is  conducted  using  the 
*RIGIDWALL_GEOMETRIC_FLAT_MOTION  option.  The  barge  bow 
mesh  is  fixed  in  space  (with  zero  displacement  boundary  conditions  speci¬ 
fied  along  the  centerline  and  the  hopper)  with  the  appropriate  approach 
angle  0  (Table  3.2)  between  the  side  of  the  barge  and  the  side  of  the 
approach  wall  (see  Figure  3.1)  for  the  run  under  study.  The  flat  nonpene¬ 
trating  rigid  wall  is  brought  into  contact  with  the  barge  at  a  constant  nor¬ 
mal  velocity  as  listed  in  Table  3.2  for  that  particular  run  until  36  in.  of 
penetration  by  the  “rigid”  (approach)  wall  into  the  barge  mesh  is 
accomplished.  The  LS-DYNA  contact  surface  option 
*CONTACT_AUTOMATIC_SINGLE_SURFACE  is  engaged  in  these 
analyses.  A  coefficient  of  friction  equal  to  0.6  (appropriate  for  steel  barge 
on  steel  armor  on  the  approach  wall)2  is  specified  for  the  contact  surface 
between  the  “rigid”  wall  and  the  barge.  The  contact  force  between  the 
approach  wall  and  the  barge  that  is  being  “crushed”  at  the  impact  corner 
region  of  the  bow  is  monitored  throughout  the  analysis.  The  nonlinear 
material  model  allows  for  a  nonlinear  structural  system  response  of  the 
multi-degree-of-freedom  structural  system  that  is  the  barge  bow. 

3.3  Results 

This  section  summarizes  the  computed  results  for  the  24  nonlinear  finite 
element  analyses.  The  deck,  the  hull  and  internal  plates,  and  the  internal 
structural  members  all  develop  elastic  and  then  plastic  strains  within  the 
steel  members  as  the  “rigid”  wall  contacts  the  barge  corner  and  then 
displaces  the  corner. 


1  Bathe  et  al.  (1999)  observe  that  crash  analyses  simulate  fast  phenomena  requiring  a  transient 
analysis  whereas  crush  analyses  simulate  slow  phenomena  that  are  approximately  analyzed  using 
static  analysis  procedures  or  implicit  dynamic  techniques.  In  recent  years,  Dr.  Gary  R.  Consolazio  and 
his  students  have  been  using  numerical  methods  to  quantify  barge  impact  loads  on  bridges  for  head- 
on  collisions  with  bridge  piers  (e.g.,  Consolazio  et  al.  2002;  Consolazio  and  Cowan  2003).  These 
calculations  have  involved  the  use  of  the  general  purpose,  nonlinear  finite  element  computer  programs 
ADINA  and  LS-DYNA.  Their  LS-DYNA  models  and  analyses  share  many  common  modeling  and  analysis 
features  with  those  discussed  in  this  report. 

2  Sullivan  (1988);  Encarta  Encyclopedia  (2004);  Grigoriev  et  al.  (1997). 
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Figure  1.4  shows  the  overhead  views  of  the  barge  impact  corner  at  initial 
contact  with  the  wall  and  after  36  in.  of  penetration  for  Case  No.  9.  Case 
No.  9  is  for  the  unusual  loading  condition  with  an  approach  angle  of  io°. 
Figure  1.5  shows  the  force  normal  to  the  approach  wall  versus  displace¬ 
ment  of  the  impact  corner  for  this  case.  This  force  represents  both  the 
normal  contact  force  applied  to  the  approach  wall  as  well  as  the  force 
imparted  by  the  wall  onto  the  multi-degree-of-freedom  structural  system 
of  the  barge.  The  normal  force  (Figure  1.5)  ranges  in  value  from  o  kips  to  a 
maximum  force  of  3,497  kips  at  4.8  in.  of  permanent  normal  deformation 
at  the  impact  corner.  The  Case  No.  9  LS-DYNA  analysis  continued  until 
there  was  36  in.  of  permanent  deformation  of  the  rigid  wall  into  the  impact 
corner  region  of  the  barge.  The  nonlinearity  in  Figure  1.5  results  from 
(1)  the  nonlinear  true-stress  versus  true-strain  relationship  for  A-36  steel 
(Figure  2.9),  and  (2)  the  removal  of  strain  integration  points  and  elements 
when  an  effective  plastic  strain  of  20%  is  achieved  within  the  shell  ele¬ 
ments  used  to  model  the  structural  members  and  plates  during  the  course 
of  the  “crushing”  analysis  of  the  elliptical  corner. 

Figure  3.2  shows  the  distribution  of  effective  plastic  strain  at  a  permanent 
deformation  of  4.8  in.  at  the  impact  corner  of  the  bow  for  Case  No.  9.  Note 
the  “crease”  or  standing  wave  of  deformation  in  the  deck  plate  in  the  upper 
right  figure  as  well  as  the  “dimpling”  of  the  elliptical  corner  at  this  level  of 
deformation.  Figure  3.3  shows  a  close-up  view  of  the  distribution  of  effec¬ 
tive  plastic  strain  at  the  impact  corner  (i.e.,  the  elliptical  corner).  Note  the 
erosion  (i.e.,  removal)  of  elements  (by  LS-DYNA)  at  the  intersection  of  the 
deck  with  the  elliptical  corner;  this  results  from  an  effective  plastic  strain 
of  0.2  (i.e.,  20%)  being  achieved  by  shell  elements  within  this  zone.  Fig¬ 
ure  3.4  shows  the  distribution  of  von  Mises  stress  at  a  permanent  defor¬ 
mation  of  4.8  in.  at  the  impact  corner  of  the  bow  for  Case  No.  9. 

Figure  3.5  shows  the  distribution  of  effective  plastic  strain  at  a  permanent 
deformation  of  36  in.  at  the  impact  corner  of  the  bow  for  Case  No.  9.  Note 
that  at  this  level  of  permanent  deformation  the  multiple  “creases”  or 
standing  waves  of  deformation  in  the  deck  plate  in  the  upper  right  figure 
as  well  as  the  “crushing”  of  the  elliptical  corner  at  36  in.  of  deformation. 
Figure  3.6  shows  a  close-up  view  of  the  distribution  of  effective  plastic 
strain  at  the  impact  corner  (i.e.,  the  elliptical  corner).  Note  the  erosion 
(i.e.,  removal)  of  elements  at  the  intersection  of  the  deck  with  the  elliptical 
corner  as  a  result  of  an  effective  plastic  strain  of  0.2  (i.e.,  20%)  being 
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Figure  3.2.  Distribution  of  effective  plastic  strain  (in  decimal  fraction)  at  4.8  in.  permanent  deformation 
and  force  versus  displacement  plot  of  the  bow  and  impact  corner  for  Case  No.  9.  Unusual  loading 

condition  with  an  approach  angle  of  10°. 

achieved  in  this  zone.  Figure  3.7  shows  the  distribution  of  von  Mises  stress 
at  a  permanent  deformation  of  36  in.  at  the  impact  corner  of  the  bow  for 
Case  No.  9.  Due  to  the  nonlinear  material  response  and  the  erosion  of 
elements,  the  normal  force  at  36  in.  of  permanent  deformation  shown  in 
Figure  1.5  is  approximately  2,300  kips,  a  significant  reduction  from  the 
maximum  force  of  3,495  kips  at  4.8  in.  of  permanent  normal  deformation 
at  the  impact  corner. 

Figures  3.8-3.13  summarize  the  force  normal  to  the  wall  versus  displace¬ 
ment  for  the  24  LS-DYNA  Cases.  The  computations  are  made  for  0-36  in. 
of  permanent  deformation  normal  to  the  wall  within  the  barge  impact 
corner  region.  The  24  cases  in  Table  3.2  are  presented  in  groups  of  four 
per  figure,  with  each  figure  sharing  a  common  load  case  (i.e.,  usual, 
unusual,  or  extreme)  and  common  approach  angle  (i.e.,  impact  angle). 
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Figure  3.3.  Distribution  of  effective  plastic  strain  (in  decimal  fraction)  at  4.8  in.  permanent  deformation 
and  erosion  of  elements  at  the  impact  corner  for  Case  No.  9.  Unusual  loading  condition  with  an  approach 

angle  of  10°. 

Figure  3.8  shows  the  force  versus  displacement  of  the  impact  corner  for 
Cases  Nos.  1-4  for  a  usual  loading  condition  with  an  approach  angle  (i.e., 
impact  angle)  of  50.  The  force  versus  displacement  results  up  to  and 
through  the  peak  range  from  3,360  to  3,605  kips.  The  peak  force  increases 
as  the  velocity  normal  to  the  wall  increases.  These  four  analyses  tend  to 
track  each  other  up  to  10  in.  of  permanent  deformation  before  there  is  a 
noticeable  divergence  in  results.  After  about  6  in.  of  permanent  deforma¬ 
tion,  there  is  a  significant  drop-off  in  force  with  additional  deformation, 
reflecting  the  nonlinear  material  response  as  well  as  the  erosion  of 
elements. 

Figure  3.9  shows  the  force  versus  displacement  of  the  impact  corner  for 
Case  Nos.  5-8  for  a  usual  loading  condition  with  an  approach  angle  (i.e., 
impact  angle)  of  10  °.  The  force  versus  displacement  results  up  to  and 
through  the  first  peak  (of  3,250-3,410  kips)  are  nearly  identical,  with  Case 
No.  5  diverging  a  little  more  from  the  others.  These  four  analyses  tend  to 
track  each  other  up  to  approximately  18  in.  of  permanent  deformation 
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Figure  3.4.  Distribution  of  von  Mises  stress  at  4.8  in.  permanent  deformation  and  force  versus  displacement 
plot  of  the  impact  corner  for  Case  No.  9.  Unusual  loading  condition  with  an  approach  angle  of  10°. 

before  there  is  a  noticeable  divergence  in  results.  After  about  6  in.  of  per¬ 
manent  deformation  there  is  a  significant  drop-off  in  force  with  additional 
deformation;  this  reflects  the  nonlinear  material  response  as  well  as  the 
erosion  of  elements.  The  authors  of  this  report  again  observe  that  there 
seems  to  be  a  dependence  of  the  nonlinear  response  on  the  (normal)  vel¬ 
ocity  specified  in  the  analysis  for  this  second  set  of  four  analyses.  Speci¬ 
fically,  Case  No.  5  seems  to  diverge  more  from  a  general  trend  at  the 
extremely  large  deformation  of  27  in.  It  cannot  be  understated  that  there  is 
significant  plastic  response  and  erosion  of  elements  at  this  stage  of  the 
nonlinear  analysis  and  the  computations  may  be  sensitive  to  small  changes 
in  plastic  responses  and  element  erosion  during  the  course  of  nonlinear 
analyses  of  the  multi-degree-of-freedom  structural  barge  bow  system.  The 
end  result  in  a  nonlinear  analysis  depends  upon  the  path  taken. 
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Figure  3.5.  Distribution  of  effective  plastic  strain  (in  decimal  fraction)  at  36  in.  permanent 
deformation  and  force  versus  displacement  plot  of  the  bow  and  impact  corner  for  Case  No.  9. 
Unusual  loading  condition  with  an  approach  angle  of  10°. 


Figure  3.10  shows  the  force  versus  displacement  of  the  impact  corner  for 
Case  Nos.  9-12  for  an  unusual  loading  condition  with  an  approach  angle 
(i.e.,  impact  angle)  of  io°.  The  force  versus  displacement  results  up  to  and 
through  the  first  peak  (of  3,450-3,500  kips)  are  identical.  The  magnitude 
of  this  first  peak  is  only  slightly  greater  than  that  computed  for  Case  Nos. 
5-8  (Table  3.3),  which  also  has  an  approach  angle  of  io°.  These  four 
analyses  tend  to  track  each  other  up  to  approximately  20  in.  of  permanent 
deformation  before  there  is  a  slight  divergence  in  results.  After  about  5  in. 
of  permanent  deformation,  there  is  a  significant  drop-off  in  force  with 
additional  deformation;  this  reflects  the  nonlinear  material  response  as 
well  as  the  erosion  of  elements.  The  authors  of  this  report  again  observe 
that  there  seems  to  be  a  dependence  of  the  nonlinear  response  on  the 
(normal)  velocity  specified  in  the  analysis  for  this  third  set  of  four  analyses 
but  with  less  divergence  than  observed  in  the  previous  eight  case  numbers. 
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Figure  3.6.  Distribution  of  effective  plastic  strain  (in  decimal  fraction)  at  36  in.  permanent  deformation 
and  erosion  of  elements  at  the  impact  corner  for  Case  No.  9.  Unusual  loading  condition 

with  an  approach  angle  of  10° 

Figure  3.11  shows  the  force  versus  displacement  of  the  impact  corner  for 
Case  Nos.  13-16  for  an  unusual  loading  condition  with  an  approach  angle 
(i.e.,  impact  angle)  of  20°.  The  force  versus  displacement  results  up  to  and 
through  the  first  peak  (of  3,180-3,390  kips)  are  identical.  The  magnitude 
of  this  first  peak  is  a  bit  lower  than  that  computed  for  Case  Nos.  9-12 
(Table  3.3),  which  has  a  smaller  approach  angle  of  io°.  These  four  analyses 
tend  to  track  each  other  up  to  approximately  25  in.  of  permanent  deforma¬ 
tion  before  there  is  a  slight  divergence  in  results,  with  Case  No.  13  showing 
more  of  a  deviation  among  the  four  analyses.  After  about  10  in.  of  per¬ 
manent  deformation  there  is  a  significant  drop-off  in  force  with  additional 
deformation;  this  reflects  the  nonlinear  material  response  as  well  as  the 
erosion  of  elements  for  this  fourth  set  of  cases.  The  authors  of  this  report 
again  observe  that  there  seems  to  be  a  dependence  of  the  nonlinear 
response  on  the  (normal)  velocity  specified  in  the  analysis  for  this  fourth 
set  of  four  analyses  and  with  about  the  same  level  of  divergence  as 
observed  in  Case  Nos.  9-12. 
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Figure  3.7.  Distribution  of  von  Mises  stress  at  36  in.  permanent  deformation  and  force  versus  displacement 
plot  of  the  impact  corner  for  Case  No.  9.  Unusual  loading  condition  with  an  approach  angle  of  10°. 


Figure  3.12  shows  the  force  versus  displacement  of  the  impact  corner  for 
Case  Nos.  17-20  for  an  extreme  loading  condition  with  an  approach  angle 
(i.e.,  impact  angle)  of  20°.  The  force  versus  displacement  results  up  to  and 
through  the  first  peak  (of  3,420-3,670  kips)  follow  the  same  general 
trend.  The  magnitude  of  this  first  peak  is  slightly  higher  than  that  com¬ 
puted  for  Case  Nos.  13-16  (Table  3.3),  which  has  the  same  approach  angle 
of  200.  These  four  analyses  tend  to  track  each  other  through  to  36  in.  of 
permanent  deformation,  with  Case  No.  20  showing  more  of  a  deviation 
among  the  four  analyses.  Past  about  11  in.  of  permanent  deformation, 
there  is  a  significant  drop-off  in  force  with  additional  deformation;  this 
reflects  the  nonlinear  material  response  as  well  as  the  erosion  of  elements 
for  this  fifth  set  of  cases.  The  authors  of  this  report  again  observe  that 
there  seems  to  be  a  dependence  of  the  nonlinear  response  on  the  (normal) 
velocity  specified  in  the  analysis  for  this  fifth  set  of  four  analyses. 
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-Case#  1 
-Case  #2 
-Case  #3 
-Case  #4 


Figure  3.8.  Force  versus  displacement  of  the  impact  corner  for  Case  Nos.  1-4. 
Usual  loading  condition  with  an  approach  angle  of  5°. 


Figure  3.13  shows  the  force  versus  displacement  of  the  impact  corner  for 
Case  Nos.  21-24  for  an  extreme  loading  condition  with  an  approach  angle 
(i.e.,  impact  angle)  of  350.  The  force  versus  displacement  results  follow  the 
same  general  trend  up  to  about  6  in.  From  6  in.  through  to  the  first  peak 
(of  3,170-3,910  kips)  there  is  a  dependency  on  the  magnitude  of  the  spe¬ 
cified  velocity  normal  to  the  approach  wall;  with  the  higher  the  normal 
velocity  resulting  in  the  higher  peak  force.  Case  Nos.  23  and  24  result  in 
the  highest  peak  forces  among  the  24  analyses  (Table  3.3). 
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- Case  #  5 

- Case  #  6 

- Case  #  7 

- Case  #  8 


Figure  3.9.  Force  versus  displacement  of  the  impact  corner  for  Case  Nos.  5-8.  Usual  loading 
condition  with  an  approach  angle  of  10°. 


Displacement,  inches 


- Case  #  9 

- Case  #10 

- Case  #  1 1 

- Case  #12 


Figure  3.10.  Force  versus  displacement  of  the  impact  corner  for  Case  Nos.  9-12. 
Unusual  loading  condition  with  an  approach  angle  of  10°. 


Table  3.3.  Summary  of  computed  maximum  force  values  for  the  24  LS-DYNA  runs. 


Category 

Impact 

Angle  0(°) 

Forward  Velocity,  Vx 

Lateral  Velocity,  Vy 

Approach  Velocity,  V 

Normal  Velocity 

Maximum  Force 

Case 

No. 

(inch/sec) 

(ft/sec) 

(inch/sec) 

(ft/sec) 

(inch/sec) 

(ft/sec) 

P(°) 

(inch/sec) 

(ft/sec) 

Force  (kips) 

Displacement 

(inch) 

Usual 

5 

6 

0.5 

0.12 

0.01 

6.0012 

0.5001 

1.1458 

0.6425 

0.0535 

3361 

2.78 

1 

Usual 

5 

6 

0.5 

1.2 

0.1 

6.1188 

0.5099 

11.3099 

1.7184 

0.1432 

3463 

3.71 

2 

Usual 

5 

24 

2 

0.12 

0.01 

24.0003 

2. 

0.2865 

2.2113 

0.1843 

3484 

3.57 

3 

Usual 

5 

24 

2 

1.2 

0.1 

24.03 

2.0025 

2.8624 

3.2872 

0.2739 

3603 

4.2 

4 

Usual 

10 

6 

0.5 

0.12 

0.01 

6.0012 

0.5001 

1.1458 

1.1601 

0.0967 

3251 

4.94 

5 

Usual 

10 

6 

0.5 

1.2 

0.1 

6.1188 

0.5099 

11.3099 

2.2237 

0.1853 

3410 

4.94 

6 

Usual 

10 

24 

2 

0.12 

0.01 

24.0003 

2. 

0.2865 

4.2857 

0.3571 

3387 

5. 

7 

Usual 

10 

24 

2 

1.2 

0.1 

24.03 

2.0025 

2.8624 

5.3493 

0.4458 

3402 

5.02 

8 

Unusual 

10 

36 

3 

4.8 

0.4 

36.3186 

3.0265 

7.5946 

10.9784 

0.9149 

3497 

4.78 

9 

Unusual 

10 

36 

3 

6 

0.5 

36.4966 

3.0414 

9.4623 

12.1602 

1.0133 

3452 

4.72 

10 

Unusual 

10 

48 

4 

4.8 

0.4 

48.2394 

4.02 

5.7106 

13.0622 

1.0885 

3474 

4.85 

11 

Unusual 

10 

48 

4 

6 

0.5 

48.3735 

4.0311 

7.125 

14.244 

1.187 

3478 

4.93 

12 

Unusual 

20 

36 

3 

4.8 

0.4 

36.3186 

3.0265 

7.5946 

16.8232 

1.4019 

3177 

9.93 

13 

Unusual 

20 

36 

3 

6 

0.5 

36.4966 

3.0414 

9.4623 

17.9509 

1.4959 

3248 

9.92 

14 

Unusual 

20 

48 

4 

4.8 

0.4 

48.2394 

4.02 

5.7106 

20.9275 

1.744 

3373 

9.8 

15 

Unusual 

20 

48 

4 

6 

0.5 

48.3735 

4.0311 

7.125 

22.0551 

1.8379 

3390 

9.87 

16 

Extreme 

20 

48 

4 

7.2 

0.6 

48.537 

4.0447 

8.5308 

23.1828 

1.9319 

3415 

9.87 

17 

Extreme 

20 

48 

4 

10.8 

0.9 

49.2 

4.1 

12.6804 

26.5656 

2.2138 

3518 

9.82 

18 

Extreme 

20 

72 

6 

7.2 

0.6 

72.3591 

6.0299 

5.7106 

31.3912 

2.6159 

3552 

9.88 

19 

Extreme 

20 

72 

6 

10.8 

0.9 

72.8055 

6.0671 

8.5308 

34.7741 

2.8978 

3674 

10.87 

20 

Extreme 

35 

48 

4 

7.2 

0.6 

48.537 

4.0447 

8.5308 

33.4296 

2.7858 

3174 

13.01 

21 

Extreme 

35 

48 

4 

10.8 

0.9 

49.2 

4.1 

12.6804 

36.3785 

3.0315 

3405 

13.82 

22 

Extreme 

35 

72 

6 

7.2 

0.6 

72.3591 

6.0299 

5.7106 

47.1954 

3.9329 

3872 

14.29 

23 

Extreme 

35 

72 

6 

10.8 

0.9 

72.8055 

6.0671 

8.5308 

50.1443 

4.1787 

3906 

14.19 

24 
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Figure  3.11.  Force  versus  displacement  of  the  impact  corner  for  Case  Nos.  13—16. 
Unusual  loading  condition  with  an  approach  angle  of  20°. 
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Figure  3.12.  Force  versus  displacement  of  the  impact  corner  for  Case  Nos.  17-20. 
Extreme  loading  condition  with  an  approach  angle  of  20°. 
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Figure  3.13.  Force  versus  displacement  of  the  impact  corner  for  Case  Nos.  21-24. 
Extreme  loading  condition  with  an  approach  angle  of  35 0 . 


Case  Nos.  21  and  22  tend  to  track  one  another,  whereas  Case  Nos.  23  and 
24  track  one  another.  All  four  analyses  tend  to  follow  the  same  general 
pattern  out  to  36  in.  of  permanent  deformation;  after  14-15  in.  of 
permanent  deformation  there  is  a  significant  drop-off  in  force  with  addi¬ 
tional  deformation,  reflecting  the  nonlinear  material  response  as  well  as 
the  erosion  of  elements  for  this  sixth  set  of  cases.  The  authors  of  this 
report  again  observe  that  there  seems  to  be  a  dependence  of  the  nonlinear 
response  on  the  (normal)  velocity  specified  in  the  analysis  for  this  sixth  set 
of  four  analyses. 
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4  Summary,  Results,  and  Conclusions 

4.1  Summary 

This  research  report  discusses  the  results  of  a  series  of  24  nonlinear  finite 
element  analysis  (using  LS-DYNA)  to  compute  the  limiting  impact  force 
due  to  the  yielding  and  buckling  of  plates  and  internal  structural  framing 
at  the  impact  corner  of  the  barge  during  its  glancing  blow  impact  with  a 
lock  approach  wall.  A  finite  element  mesh  of  the  front  half  of  the  bow  of  a 
jumbo  open-hopper  barge  consisting  of  179,238  nodes  and  177,241  ele¬ 
ments  was  constructed.  Because  of  the  localization  of  the  nonlinear 
response  during  impact  of  the  bow  of  the  barge  with  the  approach  wall, 
only  the  bow  region  of  the  jumbo  open-hopper  barge  (Figure  2.1)  is 
modeled  in  the  numerical  analyses. 

4.2  Results 

Significant  yielding  occurred  during  crushing  of  the  deck  and  hull  plates 
and  the  internal  structural  members  contained  within  the  impact  corner 
zone  of  the  bow  in  all  24  analyses.  All  24  numerical  analyses  were  carried 
out  until  36  in.  of  permanent  deformation  occurred  at  the  elliptical  corner 
region  of  the  barge.  As  such,  all  (fully  integrated)  shell  elements  in  the 
model  of  the  bow  were  specified  with  a  plastic-multilinear  material  model 
matching  the  stress-strain  data  derived  from  tests  conducted  on  A-36  plate 
steel  obtained  for  an  actual  barge.  A  large-displacement,  large-strain  form¬ 
ulation  was  used  for  the  shell  elements  in  the  LS-DYNA  impact  analyses 
and  the  true  (Cauchy)  stress,  true  (logarithmic)  strain  data  (Figure  2.7) 
specified  for  this  material  model  for  the  A-36  steel.  These  test  results  were 
obtained  from  tests  conducted  at  the  University  of  Florida  on  a  standard 
18-inch  tension  coupon  tested  by  Anderson  and  summarized  in  Consolazio 
et  al.  (2002). 

4.3  Conclusions 

This  technical  report  provides  a  basis  for  revising  the  “capping  force”  cited 
in  the  April  2004  version  of  ETL 1110-2-563.  The  capping  force  is  based  on 
the  computation  of  the  “crushing”  force  imparted  during  a  glancing  blow 
of  the  impact  corner  at  the  bow  of  the  impact  barge  with  a  lock  approach 
wall.  When  placed  in  a  severe  impact  environment  in  which  the  impact 
corner  of  the  bow  begins  to  “crush,”  the  impact  corner  at  the  bow  acts  like 
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a  structural  “fuse-plug”  and  provides  for  a  limiting  impact  force  applied  to 
the  approach  wall  by  the  barge  train.  The  results  of  these  24  nonlinear 
finite  element  analyses  of  “crushing”  of  the  barge  impact  corner  defini¬ 
tively  show  a  limiting  force  that  may  be  transferred  to  a  lock  approach  wall 
by  the  barge  train  through  the  front  barge  impact  corner. 

The  magnitude  of  this  limiting  force  ranges  in  value  from  3,170  to 
3,910  kips  for  the  24  analyses.  The  deformation  corresponding  to  this  first 
peak  force  is  dependent  upon  the  approach  angle;  it  occurred  at  per¬ 
manent  normal  deformations  of  about  5  in.  for  a  50  approach  angle;  about 
6  in.  for  a  io°  approach  angle;  about  11  in.  for  a  20°  approach  angle;  and  at 
about  15  in.  for  a  350  approach  angle. 

Table  4.1  summarizes  the  mean,  standard  deviation,  and  coefficient  of 
variation  (COV)  values  for  the  peak  forces  grouped  by  the  approach  angle 
using  the  Table  3.3  data.  Table  4.1  results  show  that  for  each  of  the  four 
approach  angle  data  groups  (i.e.,  50,  io°,  20°,  and  350),  the  average  peak 
force  values  are  within  a  narrow  range  of  3,420-3,590  kips.  The  COV  (i.e., 
the  dispersion  about  the  first  peak  mean  force  value)  is  less  than  5%  for 
approach  angles  of  20°  or  less  and  is  equal  to  10%  for  an  approach  angle  of 
35°- 


Table  4.1.  Summary  of  mean,  standard  deviation,  and  coefficient  of  variation  of  the 
computed  maximum  force  values  for  the  24  LS-DYNA  runs. 


Impact 

Angle 

Case 

First  Peak  Maximum  Force 

©C) 

No. 

Mean  Force 
(kips) 

Standard  Deviation  of 
Force  (kips) 

COV  of  Force  (kips) 

5 

1-4 

3,623 

22.470 

0.006 

10 

5-12 

3,472 

95.955 

0.028 

20 

13-20 

3,418 

161.816 

0.047 

35 

21-24 

3,589 

359.152 

0.100 
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Appendix  A:  Background— Full-Scale,  Low- 
Velocity  Controlled-lmpact  Barge 
Experiments  and  the  Development  of 
the  Empirical  Correlation 

A.l  Summary  of  the  full-scale,  low-velocity  controlled-impact 
barge  experiments 

In  December  1998,  full-scale,  low-velocity,  controlled  barge  impact  experi¬ 
ments  were  conducted  by  Patev  et  al.  (2003)  at  the  decommissioned 
Gallipolis  Lock  at  Robert  C.  Byrd  Lock  and  Dam,  Gallipolis  Ferry,  WV.  One 
of  the  many  goals  of  these  experiments  was  to  measure  the  actual  impact 
forces  normal  to  the  wall  using  a  load-measuring  device.  The  focus  of  these 
experiments  was  to  obtain  and  measure  the  baseline  response  of  an  inland 
waterway  barge,  quantify  a  multi-degree-of-freedom  system  during  the 
impact,  and  investigate  the  use  of  energy-absorbing  fenders.  The  full-scale 
experiment  used  a  297-m-long  (975-ft),  15-barge,  commercial  barge  train 
in  a  3-by-5  configuration,  as  shown  in  Figure  A.l.  Each  barge  was  a  jumbo 
open-hopper  design  10.67  by  59.45  m  (35  by  195  ft)  with  rake  barges  at  the 
front  of  the  tow.  The  barges  were  ballasted  with  anthracite  coal  to  a  draft 
of  2.74  m  (9  ft).  The  total  weight  of  the  flotilla  was  267,000  kN 
(30,012  short  tons)  with  a  total  mass  of  27,228,228  kg  (1,865.59  kips- 
sec2/ft),  which  is  equal  to  the  total  weight  divided  by  the  gravitational 
constant,  g.  A  total  of  44  impact  experiments  were  successfully  conducted 
against  the  unaltered  guide  wall  and  a  prototype  fendering  system  that  was 
attached  to  the  wall.  Approach  velocity  for  the  12  bumper  experiments 
conducted  at  the  lock,  ranged  from  268  mm/sec  to  875  mm/sec 
(0.88-2.87  ft/sec  [fps]),  with  approach  angles  (0)  ranging  from  8.8  to 
21.1 .  These  shallow  approach  angles  are  typical  of  the  approach  angles 
with  lock  approach  walls  during  glancing  blow  impacts.  For  each  of  these 
tests,  the  barge  train  was  brought  in  at  a  constant  approach  angle  0  and  at 
a  constant  velocity.  The  target  area  was  the  stiff- to-rigid  concrete  upper 
guide  wall,  lacking  the  friction-reducing  steel  armor  found  on  modern  lock 
walls.  The  approach  angle  and  velocity  for  the  12  most  credible  bumper 
experiments  are  summarized  in  Table  A.l. 
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The  load  bumper  (or,  more  specifically,  the  arc  load  beam)  used  to  record 
the  impact  force  time  histories  during  the  experiments  was  constructed  of 
mild-steel  with  an  outer  radius  of  72.6  in.,  outer  arc  length  of  43.6  in., 
cross  section  measuring  9  in.  in  width  by  5  in.  in  height,  and  separation 
between  the  6-inch  diameter  load  pins  of  35.5  in..  The  interpretation  of  the 
instrumentation  data  recorded  by  Patev  et  al.  (2003)  is  discussed  in 
Arroyo  et  al.  (2003).  The  following  summarizes  key  aspects  of  the  Arroyo 
et  al.  interpretation.  Once  the  time  of  impact  was  identified,  the  impact 
angle  (the  angle  formed  by  the  port  side  of  the  corner  barge  with  the  lock 
wall)  was  determined  from  the  global  positioning  system  (GPS)  data  cor¬ 
rected.  This  angle  is  critical  to  the  bumper  geometry  and  resulting  force 
system.  Velocity  (actually  speed)  is  simply  calculated  from  the  displace¬ 
ment  of  the  front  corner  GPS  unit  per  unit  time  (1  second).  The  initial 
orientation  of  the  bumper  relative  to  the  longitudinal  axis  of  the  barges 
was  adopted  to  be  54 0  from  the  longitudinal  axis  (local  axis  of  the  model) 
of  the  barges.  Initially,  the  recorded  forces  at  the  pins  were  assumed  to  be 
in  the  radial  direction.  The  precise  orientation  of  the  bumper  on  the  barge 
is  critical  to  this  effort.  The  as-built  orientation  of  the  bumper  was  then 
determined  from  a  combination  of  design  drawings  and  documentary 
photos.  The  survey  data  were  intended  for  this  purpose;  however,  the 
uncertainty  caused  by  the  barges  shifting  and  the  tow  drifting  against  its 
moorings  between  sightings  compromised  the  accuracy  of  these  measure¬ 
ments  sufficiently  to  make  them  unusable  for  this  purpose. 

Subsequently,  it  was  established  from  the  design  drawings  and  documen¬ 
tary  photographs  that  the  orientation  of  the  recorded  forces  was  not 
aligned  in  the  radial  direction  of  the  arc  load  beam.  Taking  into  account 
this  observed  discrepancy,  a  new  orientation  of  the  recorded  forces  was 
established.  This  second  configuration  was  analyzed  considering  the  mag¬ 
nitude  of  the  angles  associated  with  the  orientation  of  the  recorded  forces. 
The  results  of  this  analysis  indicated  that  an  impossible  geometrical 
arrangement  was  produced  by  this  second  set  of  assumptions.  A  final  geo¬ 
metrical  configuration  was  then  established  based  on:  (1)  the  range  of 
probable  angles  for  the  force  orientations  relative  to  the  radial  direction, 

(2)  the  location  of  the  bumper  related  to  the  longitudinal  axis  of  the 
barges,  and  (3)  the  appropriate  coefficient  of  friction  between  concrete  (for 
the  unarmored  wall  face)  and  (barge  corner)  steel.  Arroyo  et  al.  (2003) 
demonstrated  that  this  final  configuration  produces  reasonable  results 
based  on  the  values  of  the  coefficient  of  friction  between  the  wall  and  the 
steel  bumper  found  in  the  technical  literature,  using  the  fact  that  the 
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bumper  must  be  in  compression  during  the  impact  process.  Based  on  a 
careful  assessment  of  the  results  from  this  bumper  study,  only  eight  of  the 
initial  12  bumper  impact  experiments  were  used  in  the  empirical  correla¬ 
tion  developed  by  Arroyo  et  al.  (2003)  to  estimate  the  maximum  impact 
force  normal  to  the  wall. 

A.2  Background— empirical  correlations 

Using  values  for  the  maximum  normal  force  Fw  and  the  linear  momentum 
normal  to  the  wall,  a  best-fit  straight  line  was  calculated  for  the  eight  good 
experiments  of  1998.  This  approach  relates  the  maximum  Fw  obtained 
from  the  experiments  directly  to  the  linear  momentum.  The  least-squares 
regression  procedure  was  used  to  develop  the  best-fit  straight  line  through 
the  eight  data  points  for  the  empirical  correlation.  The  line  was  assumed  to 
start  at  the  origin  (i.e.,  no  intercept  term  was  used  for  the  linear  equation). 
The  resulting  best-fit  equation  for  this  set  of  eight  data  values  was  ( Fw)max 
-  0.435  mV  sin  0.  That  is,  a  coefficient  times  the  linear  momentum  normal 
to  the  wall  determines  the  maximum  force  normal  to  the  wall.  We  can 
observe  that  the  greater  the  magnitude  for  the  linear  momentum,  the 
larger  the  maximum  value  for  the  impact  force  normal  to  the  wall.  This 
relationship  was  based  on  low  velocity,  shallow  impact  (up  to  21. i°)  experi¬ 
ments  that,  by  definition,  do  not  account  for  factors  that  manifest  them¬ 
selves  at  higher  velocities.  Additionally,  no  damage  occurred  to  the  barge 
train  and  no  lashings  broke  during  these  eight  impact  experiments. 

The  empirical  correlation  between  the  maximum  force  normal  to  the  wall 
and  the  linear  momentum  normal  to  the  wall  immediately  before  impact 
developed  using  the  Patev  et  al.  (2003)  field  data;  development  of  empiri¬ 
cal  correlation  described  in  Arroyo  and  Ebeling  (2004)  and  Arroyo  et  al. 
(2003)  was  based  on  statistical  procedures  and  the  values  of  force 
obtained  from  the  acceptable  bumper  configuration. 

The  mathematical  form  of  Newton’s  second  law  states  that  a  resultant 
external  force  applied  to  a  body  is  equal  to  the  mass  of  the  body  multiplied 
by  the  absolute  acceleration  the  body  experiences.  Furthermore,  it  can  be 
expressed  in  terms  of  the  absolute  velocity  of  the  body  by  introducing  the 
first  derivative  with  respect  to  time  of  the  velocity,  which  is  the  accelera¬ 
tion.  One  useful  tool  that  can  be  derived  from  Newton’s  second  law, 

F  =  ma,  is  obtained  by  integrating  both  sides  of  the  equation  with  respect 
to  time.  This  integration  can  be  done  only  if  the  forces  acting  on  the  parti¬ 
cle  are  known  functions  of  time.  The  external  forces  acting  on  the  particle 
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change  the  linear  momentum.  The  mathematical  form  of  the  resulting 
expression  after  the  process  of  integration  states  that  the  impulse  during  a 
period  of  time  due  to  the  applied  impulsive  force  is  equal  to  the  difference 
in  linear  momentum  during  the  same  interval  of  time.  This  relationship 
establishes  the  Principle  of  Impulse  and  Linear  Momentum.  The  units  of 
both  impulse  and  momentum  are  force  and  time;  therefore,  impulse  and 
momentum  are  expressed  in  Newton-sec,  or  kips-sec.  The  impulsive  force 
is  a  function  of  time  and,  in  general,  varies  during  its  period  of  application. 
A  large  force  that  acts  over  a  short  period  of  time  is  called  an  impulsive 
force. 

The  linear  momentum  is  defined  as  the  mass  of  the  particle  multiplied  by 
the  velocity  of  the  particle.  It  is  a  vector  quantity  oriented  in  the  same 
direction  as  the  velocity  of  the  particle  (tangent  to  the  trajectory).  The 
velocity  of  a  barge  train  is  usually  specified  in  the  local  barge  axis:  longi¬ 
tudinal,  local  “x”  axis,  and  transverse,  local  “y”  axis.  In  this  case  two 
velocities  are  specified,  Vx  and  Vy.  To  obtain  the  velocity  normal  to  the 
wall,  an  axis  transformation  equation  is  needed.  This  expression  is  the 
following: 


K 


par 


V_ 


=  c 


-1 


Ik 


(A.l) 


where  [Cl  = 


,  and  Vpar  and  Vnorm  are  the  velocity  parallel 


cos0  -sin0 
sin0  cos0 

(global  “A”  axis)  and  normal  (global  “Y”  axis)  to  the  wall,  respectively. 
Equation  A.l  can  be  obtained  easily  from  Figure  A.2. 


Figure  A.2.  Velocity  vector  transformation,  from  local  to  global  axis  (Arroyo  and  Ebeling  2006). 
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The  empirical  correlation  between  the  maximum  force  normal  to  the  wall 
and  the  linear  momentum  normal  to  the  wall  immediately  before  impact, 
developed  by  Arroyo  et  al.  (2003)  was  based  on  statistical  procedures  and 
the  values  of  maximum  impact  force  obtained  from  the  acceptable  bumper 
configuration.  Using  values  for  the  maximum  normal  force  Fw  and  the 
linear  momentum  normal  to  the  wall,  a  best-fit  straight  line  was  calculated 
using  data  from  eight  of  the  full-scale  impact  experiments.  This  approach 
relates  the  maximum  Fw  directly  to  the  linear  momentum.  It  is  important 
to  note  that  only  one  data  point  of  the  entire  Fw  time  history  for  each  of  the 
eight  experiments  were  used  to  develop  this  empirical  correlation.  The 
least-squares  regression  procedure  was  used  to  develop  the  best-fit 
straight  line  through  the  eight  data  points  (for  the  eight  impact  experi¬ 
ments)  for  the  empirical  correlation.  The  line  was  assumed  to  start  at  the 
origin  (i.e.,  no  intercept  term  was  used  for  the  linear  equation).  The  result¬ 
ing  best-fit  straight  line,  average  -1  standard  error,  and  average  +1  stan¬ 
dard  error  lines  were  developed  and  shown  in  Figure  A.3.  The  resulting 
best-fit  equation  for  this  set  of  eight  data  values  is  ( Fw)max  =  0.435  m Vn0rm, 
with  units  of  the  resulting  force  in  kips,  mass  (including  the  mass  of  the 
loaded  barges  and  tow  boat,  but  excluding  hydrodynamic  added  mass)  in 
kips-sec2/ft,  and  approach  angle  in  degrees.  That  is,  a  coefficient  times  the 
linear  momentum  normal  to  the  wall  determines  the  maximum  force  nor¬ 
mal  to  the  wall.  We  can  observe  that  the  greater  the  magnitude  for  the 
linear  momentum,  the  larger  the  maximum  value  for  the  impact  force 
normal  to  the  wall.  This  relationship  was  based  on  low-velocity,  shallow 
impact  (up  to  21. i°)  experiments  that,  by  definition,  do  not  account  for 
factors  that  manifest  themselves  at  higher  velocities.  Additionally,  no  dam¬ 
age  occurred  to  the  flotilla  of  barges  and  no  lashings  broke  during  these 
eight  impact  experiments.  This  empirical  correlation  was  derived  using 
data  obtained  from  a  barge  train  (3-by-5)  that  had  a  velocity  normal  to  the 
wall  up  to  and  not  exceeding  173.74  mm/sec  (0.57  fps)  with  no  damage 
occurring  during  impact  events;  this  was  up  to  21.10  for  impact  angles,  and 
2,885.29  kN-sec  (between  649.84  and  1,025.48  kips-sec)  for  a  barge  train 
with  a  linear  momentum  normal  to  the  wall. 

The  maximum  normal  force  ( Fw)max  by  the  empirical  correlation  is  equal 
to  the  reaction  force  provided  by  the  lock  wall  on  the  barge  train  during 
the  impact.  Note  that  the  masses  used  to  develop  the  correlation  of  linear 
momentum  normal  to  the  wall  with  values  of  ( Fw)max  uses  the  mass  of  the 
barge  train  and  does  not  include  the  computation  of  any  hydrodynamic 
added  masses.  (However,  hydrodynamic  effects  on  the  barge  train  are 
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accounted  for  in  the  measured  impact  forces.)  A  single  lumped  mass  was 
used  to  characterize  the  barge  train  in  this  simplified  correlation. 


0  200  400  600  800  1000  1200 

Linear  Momentum  Normal  to  the  Wall  (k-s) 

(mVnorm) 


Figure  A.3.  Empirical  correlation  using  the  linear  momentum  normal  to  the  wall 
concept  (Figure  6.3  in  Arroyo  et  al.  2003). 
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Appendix  B:  Effect  of  Number  of  Processors 
on  LS-DYNA  Runs 

B.l  Cray  XT-3 

The  24  analyses  that  are  the  focus  of  this  report  were  run  on  a  Cray  XT-3 
supercomputer  located  in  the  High  Performance  Computing  (HPC)  center 
of  USACE,  Engineer  Research  and  Development  Center.  The  XT-3  is  a 
Massively  Parallel  Processor  supercomputer  that  contains  4,160  nodes, 
each  containing  one  2.6-GHz  AMD  Opteron  64-bit  dual-core  processor 
and  dedicated  memory.  Compute  nodes  comprise  4,096  of  the  nodes, 
which,  together,  have  4  GB  of  RAM.  The  remaining  64  nodes  are  service 
nodes.  The  operating  system  is  UNICOS/lc.1 

B.2  LS-DYNA  multiprocessor  capability 

LS-DYNA  is  capable  of  using  this  parallel  architecture.  When  running 
LS-DYNA,  the  user  may  specify  the  number  of  processors  to  use  in  com¬ 
pleting  the  task.  LS-DYNA  decomposes  the  model  into  a  number  of  pieces 
to  fit  the  number  of  processors  specified. 

B.3  Glancing  blow  analyses  run  times 

Analyses  are  performed  by  submitting  each  run  as  a  job  to  a  batch  queue. 
Jobs  run  as  resources  become  available.  There  are  practical  considerations 
involved  in  deciding  the  number  of  processors  to  specify  in  running  a  job. 
A  small  number  of  processors  requires  a  longer  actual  elapsed  time  for  a 
job  to  run  than  a  large  number  of  processors.  However,  the  time  a  job 
spends  waiting  in  the  job  queue  increases  as  the  number  of  processors 
requested  increases.  For  that  reason,  the  number  of  processors  used  to  run 
the  jobs  for  these  analyses  varied. 

The  24  analyses  performed  as  part  of  this  study  involved  normal  velocities 
that  spanned  a  range  of  0.6425-50.1443  in./sec.  To  obtain  a  displacement 
of  36  in.,  the  LS-DYNA  model  times  varied  from  57.6  seconds  down  to 
0.75  seconds.  A  model  time  of  0.75  seconds  completed  in  an  actual  elapsed 
time  of  31.675  hours,  using  four  processors  to  run.  Attempting  to  run  a 


1  This  information  was  obtained  from  the  HPC  Website 
http://www.erdc.hDC.mil/svstemNews/Crav  XT3/ug#configuration. 
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model  time  of  57.6  seconds  on  four  processors  would  have  required  over 
101  days.  Obviously,  that  is  not  practical.  To  obtain  results  in  a  timely 
fashion,  the  authors  decided  to  perform  analyses  using  the  number  of 
processors  shown  in  Table  B.i. 


Table  B.I.  Number  of  processors  used 
for  glancing  blow  runs. 


Case  No. 

Number  of  Processors 

1 

512 

2 

256 

3 

256 

4 

256 

5 

256 

6 

256 

7 

256 

8 

256 

9 

256 

10 

256 

11 

256 

12 

256 

13 

4 

14 

4 

15 

4 

16 

4 

17 

4 

18 

4 

19 

4 

20 

4 

21 

4 

22 

4 

23 

4 

24 

4 
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B.4  Sensitivity  of  computed  results  to  the  number 
of  processors  specified 

To  satisfy  themselves  that  LS-DYNA  results  were  not  dependent  on  the 
decomposition  of  a  problem,  the  authors  performed  analyses  of  Case 
No.  18  using  4,  32,  64,  and  256  processors.  The  computed  force  normal  to 
the  approach  wall  versus-displacement  normal  to  the  wall  results  obtained 
from  the  four  analyses  are  shown  in  Figure  B.i  for  the  fully  integrated  shell 
element  (*SECTION_SHELL;  ELFORM  16)  with  a  plastic-multilinear 
material  model  (LS-DYNA  material  No.  24, 
*MAT_PIECEWISE_LINEAR_PLASTICITY). 

Figure  B.i  shows  that  consistency  in  results  between  the  four  runs  is 
excellent  out  to  approximately  17  in.  of  permanent  deformation.  Past  that 
displacement  there  is  not  total  agreement  in  the  computed  results,  but  it 
was  concluded  by  the  authors  of  this  report  that  the  results  are  consistent, 
except  for  some  acceptable,  minor  differences.  It  cannot  be  understated 
that  there  is  significant  plastic  response  and  erosion  of  elements  at  this 
stage  of  the  nonlinear  analysis;  the  computations  may  be  sensitive  to  small 
changes  in  plastic  responses  and  element  erosion  during  the  course  of 
nonlinear  analyses  of  the  multi-degree-of-freedom  structural  barge  bow 
system.  The  end  result  in  a  nonlinear  analysis  depends  upon  the  path 
taken.  Considering  these  issues,  the  authors  are  content  that  for  the  fully 
integrated  shell  element,  the  computed  results  are  not  sensitive  to  the 
number  of  processors  specified  in  the  analysis. 

To  save  execution  time,  a  similar  parametric  study  on  the  accuracy  in 
results  to  the  number  of  processors  for  the  LS-DYNA  reduced  integrated 
shell  elements  (*SECTION_SHELL;  ELFORM  =  2;  NUMBER  OF  INTE¬ 
GRATION  POINTS  =  3)  with  hourglass  control  (^HOURGLASS;  Hourglass 
Control  Type  =  4,  Flanagan-Belytschko  stiffness  form)  was  performed. 
However,  in  plots  of  the  force  versus  permanent  displacement  normal  to 
the  wall  (similar  to  Figure  B.i),  it  is  observed  that  the  computed  results 
varied  significantly  in  the  nonlinear,  plastic  portion  of  the  contact  force 
response  if  the  number  of  processors  used  to  run  a  problem  varied.  There¬ 
fore,  only  the  fully  integrated  shell  element  model  was  specified  for  the 
computations  reported  in  this  technical  report. 
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Case  18  Number  of  Processors  Comparison 


-4  Processors 

-  32  Processors 

-  64  Processors 

-  256  Processors 


Displacement,  inches 


Figure  B.l.  Results  obtained  with  four  different  problem  decompositions  for  Case  No.  18. 
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